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The McMillan liquid crystalline models under the influence of homeo- 
tropic anchoring walls and of external fields are investigated. For 
thin systems, the existence of the critical thickness, below which the 
system does not undergo a discrete phase transition, is confirmed. 
Apparent differences between the influence of the anchoring walls 
and of external fields are elucidated by investigating the order pa- 
rameters and a temperature vs external field phase diagram for the 
bulk systems. 
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1 Introduction 

The anchoring conditions at the boundary walls strongly influence the behavior 
of the liquid crystals. The effects of the anchoring is remarkable at isotropic (I)- 
nematic(N) phase transition, since the discontinuity of extensive variables at 
I-N transition is rather small. For instance, in a thin liquid crystal system 
bounded by two parallel walls at which the liquid crystal molecules are an- 
chored homeotropically, the I-N transition occurs at higher temperature than 
the transition temperature of the bulk system 0121 El EH ■ Furthermore if the 
homeotropic anchoring is strong and the system is thin enough, the FN tran- 
sition disappears and the system undergoes a continuous change between high- 
and low-temperature phases. It is found that there is a critical thickness below 
which the system has no transition 

The effects of the anchoring at boundaries are similar to the effects of a uni- 
form external field, as indicated in the early study by Sheng pQ. Under the ex- 
ternal field which aligns liquid crystal molecules, the I-N transition temperature 
shifts higher, as under the influence of homeotropic anchoring at boundaries. 
Corresponding to the existence of the critical thickness for I-N transition in thin 
systems, there is a critical strength of the external field; i.e., the I-N transition 
vanishes and becomes a continuous change under the external field stronger 
than the critical external field ■ The similarities between the influence of 
the external field and of the anchoring boundary condition were discussed by 
Poniewierski and Sluckin 4 . In ref.@], they investigated the boundary effects 
by considering a boundary ordering potential at the walls; they compared the 
effects of external field and the effects of anchoring walls by relating the exter- 
nal field and averaged anchoring potential over the bulk. The predictions of the 
averaged anchoring potential approximation are well satisfied if the anchoring 
potential is weak. 

In the present paper, we study the influence of anchoring walls and of exter- 
nal fields upon the N and smectic A (A) phases. Then we discuss the differences 
between the influence of the anchoring conditions and of external fields. We use 
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so-called McMillan's liquid crystal model, which exhibits I, N, and A phases 0. 
We consider the situation in which the liquid crystal system is sandwiched be- 
tween two parallel walls and the liquid crystal molecules are strongly anchored 
homeotropically at the walls. Changing the temperature and the thickness of 
the system, we investigate the temperature dependence of the order parameters 
and transition behaviors. In addition we introduce effective fields due to the in- 
homogeneity of the order parameters, and compare the behaviors of the system 
under the effective fields and of the system under the external fields. 



2 Formulation 

In this paper, we use the formulation given in ref.|S] for McMillan's liquid crys- 
tal model under the external fields and under the strong anchoring boundary 
condition, which is carried out in the framework of the molecular field approxi- 
mation. 

First we introduce the self-consistent equations for a bulk system under 
uniform external fields conjugate to the order parameters. We denote the N 
and A order parameters by s and a, respectively. These order parameters 
arc defined, as was done in McMillan's original work, as s = (P2(cos0)) and 
a = (cos(27rz/d)F 2 (cos0)) where z and 6 are, respectively, the z-coordinate 
(z-axis is taken to be parallel to the director) and polar angle of a molecule. 
Here we assume that the smectic layer thickness is d. Within the mean field 
approximation, the one-body potential for McMillan model is 



V(z,cos9) = -V 



2nz 



P 2 (cost 
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where the function P 2 is the second order Legendre function P 2 (x) = (3x 2 — 1)/2, 
and the parameter a is a dimensionless interaction strength for A phase relative 
to N phase. Now we introduce two external fields h s and h a corresponding to 
the order parameters s and a, respectively. Then the self-consistent equations 
for order parameters are 



s = I(P(h a + V s),p(h a + V aa)), 
a = J{P{h s + V s),(3{h a + V aa)), 



(2a) 
(2b) 



where j3 denotes the inverse temperature. The functions I{rj, C) and J(r], C) are 
defined as 



J(77,C) = ^lnZ(ry,C), 



where Z is 

z(v,0 
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For a given set of external fields, the self-consistent equations @ have some sets 
of solutions corresponding to I, N, and A phases. Among these sets of solutions, 
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the thermodynamically stable set gives the smallest value of a function 



(3F(f3,h s ,h cr ;s,a) = 

— s + — a<J ~ ln z{oJ) ' (5) 

and the smallest value of F is the thermodynamic free energy. The other sets 
of solutions correspond to the metastable and unstable states. 

In order to treat the inhomogeneous systems, we use a discrete McMillan 
model in which the systems are divided into layers by planes parallel to the 
boundary walls. We choose the layer thickness thin enough to neglect the spa- 
tial variance of order parameters inside the layer; then we can assume the nth 
layer has its own N and A order parameters s n and a n , respectively. We assume 
that this layer thickness and the smectic layer thickness d are of the same or- 
der. In the following we restrict ourselves to the systems which have thickness 
D = (N + 2)d with integer N (Oth and N + 1th layers are in contact with 
walls). Under these assumptions, the molecular field contributions appearing 
in the right-hand side of eqs. V^s and Voaa, are divided into the intralayer 
and intcrlaycr molecular field contributions. Then the self-consistent equations 
of homogeneous systems is generalized to the self-consistent equations of inho- 
mogeneous systems as 

s„ = I{f3{h^ + V o s n ),0(h<?> + V aa n )), (6a) 
on = J(P(h™ + Vosn), (3(h?> + V aa n )), (6b) 

where Vq = Vq+2Vq'; the parameters V ' and Vq are the intralayer and interlayer 
interaction potentials, respectively. The quantities hs and hh are defined as 

= ^"(sn-i - 2s n + s n+1 ), (7a) 
= V "(<7n-i - 2a„ + a n+1 ), (7b) 

and these quantities vanish if the system is homogeneous. We shall call hs 
and h!^ the effective fields because of the similarity between (j2J and Cor- 
responding to the eq.(0, the thermodynamically stable set of solutions of eq.® 
gives the smallest value of a function 

(3F(f3,{s n },{a n }) = 

^\~ Sn + ~ aan " ln mo] [ 

Sri 

+ aV n a n a n+ i} , (8) 

n=l 



N-l 



where the sets {s n } and {<r„} are the solutions to eqs. 10. 

We impose strong homeotropic anchoring condition, i.e., the perfect nematic 
order is induced by the walls and thus so = sjv+i = 1. Furthermore we assume 
that the boundary walls are smooth planes and thus the smectic order parameter 
is also unity at each wall, i.e., ctq = &n+i = !• 
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Figure 1: Smectic A order parameter profile er„ for the systems with thick- 
ness (a) N = 12 and (b) N = 13. In both figures, curves are calculated at 
temperatures where the order parameter difference a n (T + AT) — a n (T) is max- 
imum (AT = 1CT 6 ); (a)T = 0.212192 and T = 0.212193; (b)T = 0.211898 and 
T = 0.211899. 

3 Results and Discussions 

The McMillan model exhibits qualitatively different phase behavior depending 
on the value of parameter a. The bulk system undergoes a first-order transition 
directly from I phase to A phase for a > 0.98. For a < 0.98, N phase appears 
as an intermediate phase between the I phase and A phase. Both the I-N and 
N-A transitions are first-order in 0.70 < a < 0.98. The N-A transition becomes 
a second-order transition for a < 0.70 while I-N transition remains a first-order 
transition. Although it is important to investigate systems of several a, we 
restrict ourselves here to the system which has three (I, N and A) phases and 
the two transitions between them are both first-order. We choose the parameter 
a = 0.88. The results for the other a values will be published elsewhere. 

We show in fig. ^ the order parameter profiles for systems with different 
thickness. These profiles are calculated at temperatures near the bulk transition 
temperature. Since figEfb) exhibits the first order transition and figlHa) does 
not, the critical thickness for a — 0.88 is N = 13. The transition temperatures 
of the thin system is slightly higher than the bulk transition temperature. The 
thickness dependence of the transition temperature are shown in fig [21 The 
change in transition temperature AT = T(N) — T(oo) behaves as AT ~ N~ 2 
for small N. We note that the thickness dependence of AT seems inconsistent 
with the prediction of Kelvin equation: AT ~ A -1 0E1- However, since the 
Kelvin equation is obtained by neglecting the thickness dependence of latent 
heat and of the surface tension g], the inconsistency only indicates that our 
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Figure 2: The thickness dependence of the transition temperature T(N) (cir- 
cles) for the system with parameter a — 0.88. The horizontal axis represents 
the square of inverse thickness normalized by the critical thickness N c = 13. 
The transition temperature for the bulk system is indicated with a symbol " x " . 

studies are restricted to the very thin systems where the Kelvin equation is 
not valid. Since the difference between influence of external fields and of the 
anchoring walls is mainly focused in the present study, the discussion on relation 
between our results and Kelvin equation will be published elsewhere. 

We can see the difference between the influence of external fields and of 
boundary conditions by investigating the order parameters quantitatively. The 
order parameter a n for n — 6 of thin (N = 12) system and order parameter a of 
bulk system under some external fields are shown in figCU There is a region, as 
shown in figEl where the bulk system is metastable or unstable no matter how 
the external fields are. However an order parameter of a single layer in a thin 
system pass through the region where the bulk system is not stable. In other 
words, even if we apply the uniform external fields at any strength, we cannot 
obtain the same order parameter values of a single layer. This fact clearly shows 
the difference in the influence of the external fields and of the effective fields. 

In the following, we directly compare a system under uniform external fields 
with a single layer under effective fields due to the inhomogeneity of order 
parameters, using a procedure originally introduced in ref. 9 . We focus on 
the system whose thickness is less than the critical thickness, i.e. the system 
does not have phase transition. In order for the direct comparison of external 
fields with effective fields, it is better to draw a path of a point (hi n \T), h^iT)) 
on T-h s -h a phase diagram of bulk system. However, since the three-dimensional 
phase diagram is too complicated, we consider the h^-T phase diagram, not for 
fixed h s but for h s — hs (T); i.e., we consider the phase diagram not on the 
three-dimensional T-h s -h a space but on its two-dimensional subspace defined 
by the equation h s = (T) . 

The h a -T phase diagram is shown in fig. El I n this figure we show, in addition 
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Figure 3: Dashed curves are the order parameters cr(T) for bulk systems under 
several external fields: h s = h a = 0, 2.5 x 10~ 3 , 5.0 x 10~ 3 , 7.5 x 10 -3 , and 1.0 x 
10~ 2 , from right to left. Vertical dotted lines indicate jumps of the <r(T). The 
region where the dotted lines appear cannot be reached by thermodynamically 
stable bulk systems. Solid curve is the order parameter <Jq(T) for thin (N = 12) 
system. 
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Figure 4: The temperature dependence of effective field h^\T) (n = 6 in 
the system with thickness N = 12) and the bulk phase diagram. Solid curve 
and dashed curve are the effective field and the coexisting curve between N 
and A phases, respectively. Right and left dotted curves are superheating and 
supercooling curves, respectively. 
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to the coexisting line, the superheating line and supercooling line. The N phase 
is metastable between the coexisting line and supercooling line, and the A phase 
is metastable between the coexisting line and superheating line. If a bulk system 
crosses the coexisting line of fig. 0] from high temperature to low temperature, 
the system undergoes a phase transition from N phase to A phase. The bulk 
system can change from N phase to A phase without any discrete transition 
only if the external field h a is stronger than the critical value, where the N-A 
coexisting line terminates. 

On this phase diagram, we drawn the effective field h a (T). As the temper- 
ature decreases, the h a (T) crosses the coexisting line, touches the supercooling 
line and the superheating line, and again crosses the coexisting line. During 
this process, the thin system does not undergo any discrete transition, although 
the system crosses the coexisting line. This fact shows that a simple analogy 
between the effective field and external field is not valid on discussing the phase 
transitions. Although the thin system is always in stable states, the system 
seems to have undergone metastable and unstable states during the process 
from high temperature to low temperature. This is because the stable states 
of the bulk system and a layer of the thin system are determined by different 
functions JSJ and JSJ), respectively, though these systems obey the equivalent 
equations (J2J and JJJJ, respectively. 

In summary, we have studied the influence of anchoring walls and of external 
fields on N-A phase transition, using a McMillan model. We have confirmed 
that the systems sandwiched between homeotropic anchoring walls exhibit the 
transition temperature increase and, in extremely thin systems, disappearance 
of the phase transition; these effects of anchoring walls are similar to the effects 
of external fields. However, by investigating the order parameter of a layer and 
the effective field, we have found that the influence of the anchoring walls cannot 
be understood with a simple analogy of the influence of the external field. 
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